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I. INTRODUCTION 

The precise investigation of the energy levels of hydrogenic atoms (muonium, positronium, 
hydrogen atom, deuterium, helium ions et al.) allows to obtain more exact values for 
many fundamental physical constants such as the lepton masses, the ratio of the lepton 
and proton masses, the fine structure constant, the Rydberg constant which are used for 
creating standards of units 0. The insertion of new simple atomic systems in the range 
of experimental investigation can lead to significant progress in solving of these problems. 
The measurement of muonic hydrogen Lamb shift at PSI (Paul Sherrer Institute) with a 
precision of 30 ppm will allow to improve our knowledge of the proton charge radius by an 
order of the magnitude |2] . Another important problem is connected with the measurement 
of the ground state hyperfine splitting (HFS) in muonic hydrogen 0, 0]- In the case of 
electronic hydrogen HFS was measured with extremely high accuracy many years ago 

AE^^s = 1420 405 751.7667(9) kHz. (1) 

The corresponding theoretical expression of the hydrogen hyperfine splitting can be writ- 
ten in the form {AE^hfs = MAuf^^g) jQ]: 

^E'h'fS = Ml + 5^^^ + + E, = 8 4 /^P^P^e (2) 
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where /ip is the proton magnetic moment, me, rUp are the masses of the electron and proton. 
The calculation of different corrections to the Fermi energy Ep has a long history. Modern 
status in the theory of hydrogenic atoms was presented in details in S^^^ denotes the 
contribution of higher-order quantumelectrodynamical effects. Corrections 6^ and 6^ take 
into account the influence of strong interaction. 6^ describes the effects of proton finite- 
size and recoil contribution. 6^ is the correction due to the proton polarizability. Basic 
uncertainties of theoretical result (2) are related with 5^ and 6^. 
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FIG. 1: Two-photon Feynman amplitudes determining the correction of the proton polarizability 
in the hyperfine splitting of the hydrogen atom. 

The main part of the one-loop proton structure correction is determined by the following 
expression (the Zemach correction) |7]: 



7r2 J p4 



Ge(p^)G'm(p 
1 + K 



1 



EFi-2fia)Rp. (3) 



The Zemach radius Rp is determined by the densities of electric charge and magnetic mo- 
ment. It is considered as a fundamental parameter of the proton structure along with the 
proton ch arg e radius. There exist three possibilities to determine the numerical value Rp 
|1, [13, El- One approach based on the analysis of the world data on e — p scattering 
gives Rp = 1.086 ± 0.012 fm j^. Another method uses the comparison of theoretical and 
experimental results (l)-(2) for the hydrogen atom. In this case the value -Rp = 1.043(16) fm 
is obtained in j^. The third method can be based on the comparison of the future experi- 
mental data and theoretical result for muonic hydrogen 0, 0, . The proton polarizability 
correction is important among other contributions of order 10~^. Numerical estimation of 
6^ obtained in 1^ serves at present as a reliable guide for defining the total value of the 
HFS in the hydrogen BSlHini- 

The aim of this work consists in the investigation of the proton polarizability correction in 
the hydrogen HFS. We performed new calculation of 6^ using the isobar model describing 
the processes of photo- and electroproduction of vr, rj mesons, nucleon resonances on the 
nucleon in the resonance region, and on evolution equations for the parton distributions in 
deep inelastic region. 
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II. GENERAL FORMALISM 



The main contribution to 5^ is determined by two-photon diagrams, shown in Fig. 1. 
The corresponding amphtudes of virtual Compton scattering on the proton can be expressed 
through nucleon polarized structure functions Gi(z/, Q^) and G'}.{ p,Q ^). Inelastic contribu- 
tion of the diagrams (a), (b) Fig. 1 can be presented in the form [l^L Ilii llfil \v\ : 



where i/^/i determines the pion-nucleon threshold: 



z^t/i = H , (7) 



2mp 



and the functions f3i 2 have the form: 



(3,(6) = 36 - 26' -2(2-6) J 6(6 + 1), (8) 



= 1 + 26- 2^6(6 + 1), 6 = v'IQ\ (9) 

-^2(Q^) is the Pauli form factor of the proton, k is the proton anomalous magnetic moment: 
K=l. 792847351(28) 

The polarized structure functions g\(v, Q'^) and (72 Q"^) enter in the antiswnmetric part 
of hadronic tensor W^,^, describing lepton- nucleon deep inelastic scattering |lq : 

w,. = W'-i? + wl'iK (10) 



(12) 



where e^^^a/s is the totally antisymmetric tensor in four dimensions, Q^) = 

mpZ/G'i(z/, Q^), g2(^,Q'^) = fnp^'^G2(^,Q'^), P is the four-momentum of the nucleon, x = 
Q^/2mpiy is the Bjorken variable, S is the proton spin four-vector, normalized to S'^ = — 1, 
= —Q^ is the square of the four-momentum transfer. The invariant quantity P ■ g is 
related to the energy transfer u in the proton rest frame: P-q = TUpV. The invariant mass of 
the electroproduced hadronic system, W, is then W'^ = m'^ + 2mpU — Q'^ = m'^ + Q'^ (1 / x — 1). 
Here Wi and W2 are the structure functions for unpolarized scattering. In the DIS regime 
the invariant mass W must be greater than any resonance in the nucleon. The threshold 
between the resonance region and the deep-inelastic region is not well defined, but it is 
usually taken to be at about W'^ = 4 GeV'^. 
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Hadronic tensor Wf^^y is proportional to the imaginary part of the off-shell Compton 
amplitude for the forward scattering of virtual photons on nucleons: 7*A^ 'y*N. The 
photon-nucleon interaction depends on the photon polarization as well as on the nucleon 
one. This gives four independent helicity amplitudes of the form Mat-cd, with a, b, c, d 
values for the helicities of the photon and nucleon initial and final states: 



^l,l/2;l,l/2? ^1 -l/2;l,-l/2, Afo,l/2;0,l/2, Afi^i/2;0,-l/2- 



These components correspond to the four structure functions Wi, W2,gi,g2- All other pos- 
sible combinations of initial and final photon and nucleon helicities are related to the above 
by time reversal and parity transformation. 

The proton spin structure functions can be measured in the inelastic scattering of polar- 
ized electrons on polarized protons. Recent improvements in polarized lepton beams and 
nucleon targets have made it possible to make accurate measurements of nucleon polarized 
structure functions ^1,2 in experiments at SLAG, CERN and DESY HI H H 13, H . 
The spin dependent structure functions can be expressed in terms of virtual photon- 
absorption cross sections E~ 



1712 ■ K 



87r2a(l + QVz/2 



cri/2(i^, Q^) - 0-3/2(1/, Q^) + 



(13) 



7712 ■ K 



2z/ 



-0-1/2(1/, Q ) + CT3/2(Z/, Q ) + ^=CrTL(z/, Q 



(14) 



where K is the kinematical flux factor for virtual photons, (T1/2, (^3/2 are the virtual pho- 
toabsorption transverse cross sections for the total photon-nucleon helicity of 1/2 and 3/2 
respectively, cttl is the interference term between the transverse and longitudinal photon- 
nucleon amplitudes. In this work we calculate contribution AEf^pg on the basis of the latest 
experimental data on the structure functions 5'i,2(^, Q^) and theoretical predictions for the 
cross sections <Ji/2,3/2,tl- 

The proton polarizability contribution to HPS in the resonance region is determined by 
the processes of photo- and electroproduction on nucleons of the pions and some prominent 
baryon resonances. The amplitudes of such reactions are shown in Fig. 2. 

To obtain correction (4) in the resonance region {W"^ < AGeV^) we use the Breit-Wigner 
parameterization for the photoabsorption cross sections in Eqs.(13)-(14), suggested in (25l.l26l 
|2^, ^ ^ ^ . In the considered region of the variables k"^ , W the most contribution is 
given by five resonances: P33(1232), 5ii(1535), ^13(1520), Pii(1440), Fi5(1680). Accounting 
the resonance decays to the Ntt— and Nrj— states we can express the absorption cross 
sections ai/2 and (J3/2 as follows: 



'7l/2,3/2 




■Nn 



Arrir 



{W^-Miy + W^Tl,MnTn 



/2,3/2| 



(15) 



where Ai/2,3/2 are transverse electromagnetic helicity amplitudes. 




X = 0.3 GeV. 



(16) 
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FIG. 2: Feynman amplitudes for the proton polar izability correction in the resonance region. 
Solid, double solid, wave and dashed lines correspond to nucleon, baryon resonance, photon and 
pion correspondingly. 

The resonance parameters Tr, Mji, ji, j2, Tfot are taken from jlZ) H, IE 23|- In accordance 
with 27, 2^, 34 1 the parameterization of one-pion decay width is 




Tn^NAq) = ( - ) ( ) ,C = 0.3 GeV (17) 

for the P33(1232) and 




fe'-^H^^j [f^) , (18) 

for Di3(1520), Pii(1440), Fi5(1680). / is the pion angular momentum and 5^ = (M^j— 
rrip — rriTr)'^ + r|j/4. Here q {k) and qji {kji) denote the c.m.s. pion (photon) momenta of 
resonances with the mass M and Mr respectively. In the case of 511(1535) we take into 
account tiN and rfN decay modes (29l . l34| : 

2 

I R^TT,,? - — OTT.r,! R o ,^2 ' y^^) 

where b^^ j^ is the tt [rj) branching ratio. 

The cross section axL is determined by an expression similar to Eq.(15), containing the 
product (5*^/2 ■ ^1/2 + ^1/2^1/^) [19] • The calculation of helicity amplitudes A1/2, ^3/2 
and longitudinal amplitude St/'?., as functions of Q^, was done on the basis of constituent 
quark model (CQM) in ji^, l36l 1371 l38L l39l. li^ . In the real photon limit = we take 
corresponding resonance amplitudes from j33|. 

The two-pion decay modes of the higher nucleon resonances (5*11(1535), Z^i3(1520), 
Pn(1440) and Fi5(1680) were described phenomeno logically using two-step process as in 
[27[. The high-lying nucleon resonance R can decay first into A^* (P33(1232) or Pii(1440)) 
and a pion or into a nucleon and p- or a-meson. Then the new resonances decay into a 
nucleon and a pion or two pions: 



R ^ r + a 




iV + TT + TT, 

+ TT + vr. 
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FIG. 3: The proton polarized structure function Q^) as the function of variables (0 1) 

Gev'^ and W (1.1 ^ 2.0) GeV . 

The total decay width of such processes can be presented as a phase space weight integral 
over the mass distribution of the intermediate resonance r = N*, p, a {a = n, N): 

R^r+a{VV } - ^ (^n^ ^^2 _ ^2)2 + ;,2r2^^^(^) (1^ - ^3 - 2m^)2 + ' ^^^^ 

where C = 0.3 GeV^, the factor must be taken from the constraint condition: 
r R_^j.+a{WR) coincides with the experimental data, P/ is the three momentum of the reso- 
nance r in the rest frame of R. Tr^tot is the total width of the resonance r. The decay width 
of the meson resonance is parameterized similarly to that of the P33(1232): 

(\ 2Jr+l 2 I x2 



where and /i are the mean mass and the actual mass of the meson resonance, q and 
are the pion three momenta in the rest frame of the resonance with masses /i and m^, Jr 
and Tr are the spin and decay width of the resonance with the mass rrir. 

Main nonresonant contribution to the cross sections ct^.l in the resonance region is de- 
termined by the Born terms constructed on the basis of Lagrangians of jMN, 'jittt, nNN 
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FIG. 4: The proton polarized structure function g2{W,Q'^) as the function of variables (0 ■ 
Gev'^ and W (1.1 ^2.0) GeV. 



1) 



interactions. Another part of nonresonant background comprises the t- channel contribu- 
tions of p, uj mesons obtained by means of effective Lagrangians 77?^^, VNN interactions 
{y = p,uj) [sj. In the unitary isobar model accounting the Born terms, the vector meson, 
nucleon resonance contributions and the interference terms we calculated the cross sections 
c"t,l by means of numerical program MAID (http://www.kph-uni-mainz.de/MAID) in the 
resonance region as the functions of two variables W and Q^. The obtained nucleon polar- 
ized structure functions (71,2 (M^, Q^) are presented in Figs. 3-4. These results for the structure 
function gi{W, Q"^) are in qualitative agreement with experimental data. The particular sig- 
nificance in the study of the spin- dependent properties of baryon resonances belongs to 



Gerasimov-Drell-Hearn (GDH) sum rule |41 



Ami 



00 djy 

— [0-1/2 (z^,0) - 0-3/2(2/, 0)]. 
i^th ^ 



(22) 



The GDH sum rule rests on the basic physical principles and an unsubtracted dispersion 
relation applied to the forward Compton amplitude. 

To construct the spin-dependent structure functions of the proton in deep inelastic region 
we can use the evolution equations for the quark and gluon distributions 4^ |: 



dqi{x,Q'^ 



2^ 



y 



X 



(23) 
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dg{x,Q^) 



2-K 



y 



Y.q.{y,Q')P,,[-^]+g{y,Q')P,, 



(24) 



where the sum is considered over all quarks and antiquarks. Pqq, Pgg, Pgg, Pgg are the 
quark-gluon splitting functions [i^ . Numerical solution of the integrodifferential evolution 
equations (23), (24) by means of the method suggested in ji^] allows to obtain the parton 
distributions and the structure functions (71,2(2;, Q"^) for different values of a photon momen- 
tum squared Q^. Corres p ondin g numerical results are in good agreement with the world 
experimental data 



Uorresp 
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III. NUMERICAL RESULTS 

In this paper we calculate the proton polarizability correction to the hyperfine splitting 
of the ground state in the hydrogen atom on the basis of the isobar model describing the 
processes of low-energy scattering of virtual photons on nucleons and the evolution equations 
for the parton distributions. These two significant ingredients of the calculation allow to 
construct the absorption cross sections for transversely and longitudinally polarized photons 
by nucleons and to express the structure functions gi,2iW,Q'^) (13), (14), which determine 
required contribution (4). 

The values of contributions 6[, 62 and the total contribution 6^, obtained after the 
numerical integration in the resonance and nonresonant regions are as follows in electronic 
and muonic hydrogen correspondingly: 

5f = 2.6 ppm, 62 = -0.4 ppm, 5^ = (2.2 ± 0.8) ppm, (25) 

(5f = 5.18 X 10-\ (5f = -0.48 x 10~^ 6^ = (4.70 ± 1.04) x 10"^ (26) 

The difference of obtained number for electronic hydrogen from the result of our previous 
work (1.4 0.6) ppm is connected with new contributions (nonresonant background in 
the resonance region and two pion decays of the resonances) which we considered in this 
study. There exists a number of theoretical uncertainties associated with quantities entering 



in the correction (4). In the improved isobar model [3l|, |45|, |46| containing 14 resonances, 
we can omit theoretical error which arises due to the insertion of other high-lying nucleon 
resonances. On our sight the main theoretical error is closely related with the calculation 
of the helicity amplitudes Ai/2{Q'^), A3/2{Q^), Si/2{Q'^) in the quark model based on the 
oscillator potential [18]. Only systematical experimental data for the helicity amplitudes 
of the photoproduction on the nucleons y4i/2(0), ^3/2(0) are known with sufficiently high 
accuracy to the present js^. In the case of amplitudes for the electroproduction of the 
nucleon resonances experimental data contain only their values at several points Q^. So, 
we have no consistent check for the predictions of the oscillator model. Possible theoretical 
uncertainty connected with the calculation of amplitudes yli/2(<5^), ^3/2 (Q^), Si/2{Q'^) with 
the account of relativistic corrections can attain the value of order 10%. Then the theoretical 
error for the correction (4) in the resonance region comprises 20% from the obtained value. 
We solved DGLAP equations in the NLO approximation, so possible uncertainty in 5^ can 
comprise near 10% of obtained result in the nonresonant region. The other source of the 
theoretical uncertainty arises from the experimental data errors in the <1 GeV^ region. 
We estimated it at a level of about 20% of the contribution 6^ at < 1 GeV'^ in the 
nonresonant region. 



9 




1100 1200 1300 1400 1500 

W (MeV) 




-2.5 - 1 

1100 1200 1300 1400 1500 1600 

W (MeV) 

F i 




1100 1200 1300 1400 1500 1600 

W (MeV) 



FIG. 5: Plots of 7r+ virtual photoproduction cross sections at = 0.4 GeV as the functions of 
the W (1100-h 1600) MeV, obtained on CLAS detector. The solid curves correspond to the results 
from MAID. 
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New experimental data for the electroproduction cross sections in the reaction ep e'vr+n 
in the resonance region were obtained recently on CLAS detector ji^l • Their comparison with 
the calculations carried out on the basis of the unitary isobar model (MAID) is presented in 
Fig. 5. It evidently shows that MAID, which we used in our calculation, gives the numerical 
values for the electroproduction cross sections in the regions of the A isobar and resonances 
Dii(1520), 5*11(1535), which are slightly higher (approximately by 5%) than the experimental 
data. 

The work was supported by the Russian Foundation for Basic Research (grant No. 06- 
02-16821). 
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